Lattice analysis for the energy scale of QCD phenomena 
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We formulate a new framework in lattice QCD to study the relevant energy scale of QCD phe- 
nomena. By considering the Fourier transformation of link variable, we can investigate the intrinsic 
energy scale of a physical quantity nonperturbatively. This framework is broadly available for all 
lattice QCD calculations. We apply this framework for the quark-antiquark potential and meson 
masses in quenched lattice QCD. The gluonic energy scale relevant for the confinement is found to 
be less than 1 GeV in the Landau or Coulomb gauge. 
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Physics covers a huge number of phenomena on vari- 
ous scales, i.e., from elementary particles to the universe. 
Nature is described in different ways depending on the 
energy scale. The energy scale characterizes physical phe- 
nomena and is important in understanding them. 

The classical QCD Lagrangian is scale invariant, ex- 
cept for quark mass terms. However, quantum effects 
yield the scale anomaly and violate the scale invariance, 
and the energy scale appears in the real world. The ap- 
pearance of the energy scale enriches QCD phenomena 
and creates their hierarchy. In high energy, due to the 
weak coupling or asymptotic freedom, the phenomena 
are governed by perturbative QCD [3, 0|- In contrast, 
in low energy, there are complicated and important non- 
perturbative phenomena, i.e., quark confinement, chiral 
symmetry breaking, and so on [!, 0, |{|. Many of the 
QCD phenomena would have their intrinsic energy scale. 
These phenomena originate from a single Lagrangian but 
seem like different physics. The knowledge about the en- 
ergy scale of QCD phenomena would be useful in under- 
standing complicated QCD dynamics. 

One possible way to study the intrinsic energy scale 
is to calculate a physical quantity after cutting a part of 
the region of momentum space artificially. From the rela- 
tion between the cut region and the resulting quantity, we 
can understand the relevant energy scale for the quantity. 
This concept is used, for example, in the Swinger-Dyson 
approach [g]. In this letter, we introduce this concept 
to lattice QCD, which is the nonperturbative and first- 
principle calculation 0, @|. The fundamental degree of 
freedom in lattice QCD is the link variable, which rep- 
resents the gluon on lattice. We cut a part of the link 
variable in momentum space, and investigate its effect on 
the resulting quantity. 

The procedure is formulated as the following five steps. 

1. We generate a gauge configuration, i.e., link vari- 
ables C/ m (x)(g SU(3)), by the usual Monte Carlo simu- 
lation in lattice QCD. Since the procedure is not gauge 
independent, we fix the configurations with some gauge. 
It is desirable to choose a physically meaningful gauge 
on lattice, such as the Landau gauge, which has the cor- 
responding continuum theory and the minimum fluctua- 



tion of the gauge field. In this letter, we mainly use the 
Landau gauge for numerical calculation, unless otherwise 
stated. 

2. Using discrete Fourier transformation in the 4- 
dimensional Euclid space, we construct the momentum- 
space link variable U^p) as 

= £i exp^X^av), (1) 

X 

where L is a number of lattice sites in one direction. 
When numerical simulation is performed in a certain 
lattice size with periodic boundary condition, the cor- 
responding momentum space is the lattice with the same 
structure. The first Brillouin zone of the momentum 
space is a hypercube with each side of (— iv/a, it /a]. 

3. We impose a "cut" on U fl (p) in some region of the 
momentum space, for example, a ultraviolet (UV) cut or 
a infrared (IR) cut. Outside the cut, the link variable is 
replaced with the free link variable UJ C0 (p) = 6po- In this 
letter, we estimate the energy scale by the 4-momentum 
length \fp^ = ^J^^P^Pfj., and define the momentum- 
space link variable Uj}{p) with the UV cut Atjv and the 
IR cut Air, as 

{ °po WP > A uv or VP < Air) 

These cuts are schematically depicted in Fig. [1] Of 
course, we can consider other kinds of cuts, for exam- 
ple, the cut for the intermediate-momentum region, the 
anisotropic cut, and so on. 

4. We reconstruct the coordinate-space link variable 
U^(x) with the cut by the inverse Fourier transforma- 
tion. Since the naive inverse Fourier transformation of 
U£(p) is not a SU(3) matrix, we define U£(x)(e SU(3)) 
by maximizing 

Rm[Ut(x)%(x)], (3) 

where 

U » = E ^» expHE^av)- (4) 
p 
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FIG. 1: The two-dimensional schematic figure of the UV cut 
Auv and the IR cut Air on momentum-space lattice. The 
momentum-space link variable Up(p) is replaced with the free 
one Ujf ee (p) in the cut (shaded) regions. a p is the momentum- 
space lattice spacing. 



5. Using this link variable U^(x), we compute physical 
quantities in just the same way as the original SU(3) 
lattice QCD. 

By this procedure, we can investigate what energy 
scale of the gluon is essential for a QCD quantity non- 
perturbatively and quantitatively. We have only to re- 
place Uj}(x) instead of U^(x) in the usual lattice QCD 
calculation. Therefore this framework can be applied to 
whatever can be calculated by lattice QCD, and would 
be widely useful for QCD phenomena. 

We would like to apply this procedure to several fun- 
damental quantities in quenched lattice QCD. The nu- 
merical simulation is performed with the f3 — 6.0 plaque- 
tte gauge action on 16 4 lattice. The lattice spacing a is 
about 0.10 fm 0]. Then the corresponding momentum- 
space lattice is also 16 4 lattice, and the momentum-space 
lattice spacing a p — 2n/La is about 0.77 GeV. 

First, we calculate the quark-antiquark (QQ) potential 
with the UV or IR cut. The QQ potential is extracted 
from the expectation value of the Wilson loop. To en- 
hance the ground-state component, the smearing method 
is used for the spatial link variables in the Wilson loop 
[l0( . The configuration number is 50, but effectively large 
by averaging all the parallel-translated Wilson loops in 
one configuration. 

The physical QQ potential is written with a Coulomb 
plus linear potential as 



A 

V{R) =crR--+C, 
R 



(5) 
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FIG. 2: The QQ potential V(R) with the IR cut Am against 
the interquark distance R. The unit is scaled with the lattice 
spacings a ~ 0.10 fm and a p ~ 0.77 GeV. The broken line is 
the original QQ potential in lattice QCD. 



pcrturbativc one-gluon exchange, and its coefficient A is 
about 0.26. The linear confinement potential is purely 
a nonperturbative phenomenon. Its coefficient, or the 
string tension, a is about 0.89 GeV/fm, and aa 2 ~ 0.051 
in = 6.0 lattice unit. The constant term C is physically 
irrelevant and depends on the rcgularization. 

We show the QQ potential V(R) with the IR cut Air 
in Fig. [2] When AiR/a p = 1 and 1.1, the string ten- 
sion slightly decreases from the original QQ potential. 
When Ajr/ a p > 1.5, the linear conffiicmcnt potential 
completely disappears and the QQ potential becomes 
only the perturbative Coulomb potential. In Table [H 
we list the asymptotic string tension a asym which is esti- 
mated by fitting the QQ potential in 4 < R/a < 8 with 
a linear function a aS y m R + const. The asymptotic string 
tension is approximately zero in Am/dp > 1.5. This re- 
sult means that the QQ confinement potential is made 
from the low-momentum gluon, and that its physical en- 
ergy scale is less than 1 GeV in the Landau gauge. 

Since the procedure is not gauge invariant, it is desired 
to check the gauge dependence of the resulting energy 



TABLE I: The asymptotic string tension <T asym with the IR 
cut Air. The results of the Landau gauge and the Coulomb 
gauge are shown. In original lattice QCD with f3 = 6.0, the 
string tension aa 2 is about 0.051. 



Am/op 



(Landau) 



(Coulomb) 



where R is the distance between the quark and the an- 
tiquark llj. The Coulomb potential is explained by the 



1.0 
1.1 
1.5 
2.0 
3.0 



0.0469(58) 
0.0311(49) 
-0.0019(20) 
-0.0132(6) 
0.0003(12) 



0.0289(58) 
0.0190(59) 
0.0024(25) 
-0.0041(10) 
0.0058(10) 
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scale among several reasonable gauge choices in the step 
1 . In Table HI we also show the result of the Coulomb 
gauge for the gauge choice. The value of <7 asym is differ- 
ent between the Coulomb gauge and the Landau gauge. 
However, the qualitative behavior is the same, and the 
energy scale for the confinement potential is less than 
about I GeV also in the Coulomb gauge. 

The QQ potential with the UV cut Auv is shown in 
Fig. [3] As the UV cut becomes lower, the Coulomb po- 
tential gradually weakens. In contrast, the string tension 
remains unchanged regardless of the UV cut. As for the 
constant term C in Eq. ([5]), it is considered to be mainly 
given by the lattice regularization of the Coulomb singu- 
larity, that is, a perturbative origin. Actually, the con- 
stant term becomes smaller with lower UV cut. With 
the very low UV cut, the QQ potential is almost the lin- 
ear potential, and the Coulomb coefficient and the con- 
stant term is nearly zero or slightly negative. For exam- 
ple, with Ayv/flp = 2, the fitting result by Eq. (J5]) is 
era 2 = 0.0459(15), A = -0.056(5), and Ca = -0.095(7). 

Now we consider another type of the cut, i.e., the cut 
by max(pi,p2)P3)P4) instead of \fp^ in Eq. ([2]). This 
type of the cut forms a 4-dimensional hypercube, and 
respects the momentum-space lattice structure. We find 
that the QQ potential with this type of the cut displays 
the same qualitative behavior as before in both the UV 
and IR case, and gives almost the same energy scale. A 
part of the results is shown in Fig. [U 

Next, we apply this framework to the meson mass cal- 
culation. The quark propagator is calculated with the 
clover fermion action, which is the 0(a)-improved Wilson 
fcrmion action. The hopping parameters are k = 0.1200, 
0.1300, and 0.1340, and the corresponding pion masses 
are about 2.9 GeV, 1.8 GeV, and 1.3 GeV, respectively 
Q. The configuration number is 100 here. 

Figure [5] shows pion mass m T and p-meson mass m p 
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FIG. 3: The QQ potential with the UV cut Auv- The nota- 
tion is the same as Fig. [2] 
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FIG. 4: The QQ potential V(R) with the IR cut Air by 
max(pi,p2,P3,P4) instead of \f[P in Eq. ([2]). The notation is 



the same as Fig. [2] 



with the IR cut Air. When more IR region is cut, the 
pion mass and the p-meson mass become smaller. In 
addition, the pion and p-meson masses degenerate in 
AiR,/a p > 2. This energy scale is consistent with the 
energy scale of the disappearance of the QQ confinement 
potential. Therefore this degeneracy is considered to sug- 
gest that the quark and the antiquark become unbound 
and quasi-free and that m v and m p are equal to the sum 
of the quark and antiquark masses. We have also calcu- 
lated these meson masses with the UV cut, and found 
that these masses decrease also by the UV cut but the 
degeneracy does not occur. It would be also interesting 
to investigate the meson masses and the chiral conden- 
sate by this procedure with the staggered fermion or the 
overlap fermion. 

In summary, we have proposed the new framework in 
lattice QCD to investigate the relation between a QCD 
phenomenon and its typical energy scale. Using Fourier 
transformation, we construct momentum-space link vari- 
ables, cut a part of them in momentum space, and re- 
construct coordinate-space link variables with the cut. 
By lattice QCD calculations with this link variables, we 
can investigate the gluonic energy scale of QCD quanti- 
ties. Since what is needed is only the replacement of the 
link variables, this is easily available for all lattice QCD 
calculations. In addition, this kind of framework would 
be broadly applied to not only lattice QCD but also any 
discrctized theory By such a framework, we can sys- 
tematically analyze the momentum structure of physical 
phenomena. 

We have applied this framework to the calculation for 
the QQ potential and meson masses. When IR gluons are 
cut, the nonperturbative confinement potential vanishes 
in the QQ potential, and pion and p-meson degenerate 
and seem to be the quasi-free quark and antiquark pair. 
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A IR /a p 

FIG. 5: Pion mass rthr and p-meson mass m p with the IR cut 
Air. The results of three hopping parameters k = 0.1200, 
0.1300, and 0.1340 are shown. The unit is scaled with the 
lattice spacings a ~ 0.10 fm and a p ~ 0.77 GeV. The values 
at Am/a p = is calculated in original lattice QCD 



The relevant gluonic energy scale for quark confinement 
is less than about 1 GeV in the Landau gauge and the 
Coulomb gauge. 
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